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THE  WAVE  RESISTANCE  OP  BODIES  OP  REVOLUTION 


by 

Georg  P.  Welnblum,  D.Eng. 

ABSTRACT 

Following  a  brief  review  of  prior  work  on  wave  resistance  of  bodies  of  revo¬ 
lution  carried  out  by  Havelock  and  Weinblum  a  discussion  is  presented  of  the  appro¬ 
ximate  relations  between  the  shape  of  sectional-area  curves  and  of  hydrodynamic 
irregularity  distributions.  The  latter  are  expressed  by  polynomials,  which  lend  them¬ 
selves  to  an  evaluation  of  the  basic  resistance  integrals  by  computing  Intermediate 
integrals.  Values  at  the  functions  thus  obtained  are  tabulated  in  an  appendix.  These 
functions  are  then  used  to  calculate  the  resistance  of  some  simple  bodies  of  revolu¬ 
tion.  Also  investigated  is  how  the  resistance  is  influenced  by  asymmetry  with  respect 
to  midship  section.  Distributions  leading  to  bodies  of  least  wave  resistance  are  cal¬ 
culated,  assuming  rather  severe  restrictions.  A  rather  complete  set  of  resistance 
curves  is  given  for  an  important  family  of  bodies. 


1 .  INTRODUCTION 

When  a  body  moves  uniformly  and  rectilinearly  in  an  unbounded 
liquid  the  only  resistance  experienced  by  it  is  the  viscous  drag.  Our 
knowledge  as  to  how  this  drag  depends  upon  the  body  form  is  very  limited, 
but  it  is  well-established  that  for  streamlined,  elongated  hulls— with  which 
we  are  only  concerned — the  drag  is  roughly  proportional  to  the  wetted  sur¬ 
face  and  is  rather  insensitive  to  reasonable  changes  in  the  shape.1*  The 
well-known  airship  form  with  a  rather  blunt  forebody  and  finer  tail  appears 
to  be  close  to  the  minimum  resistance  attainable,  although  it  must  be  empha¬ 
sized  that  earlier  resistance  data  obtained  in  wind  tunnels  at  low  Reynolds 
numbers  are  utterly  unreliable.  But  that  there  is  a  slight  advantage  in 
introducing  some  asymmetry  with  respect  to  the  midship  section  appears  to 
be  unquestioned,  at  least  when  larger  end-radii  are  used.  Matters  become 
different  when  a  body  moves  close  to  the  free  surface;  see  Figure  1 .  A 
wave  pattern  is  then  produced  and  therefore  a  wave  resistance  arises.  The 
laws  governing  the  wave  resistance  R  are  quite  different  from  those  valid 


1References  are  listed  on  page  $8. 

♦Problems  of  cavitation  are  not  considered  here. 


for  the  viscous  drag  Rv«  Hence,  in 
this  case  forms  of  least  total  resis¬ 
tance  R.  must  be  derived  from  addi- 

f  z 

tional  considerations  and  may  differ, 
at  least  in  principle,  from  the  fa¬ 
miliar  streamlined  forms. 

In  the  present  report  it 

Pigure  1  -  Scheme  of  Submerged  Body  intended  to  analyze  the  wave  re¬ 

sistance  of  a  rather  wide  class  of 

elongated  bodies  of  revolution,  using  an  integral  relation  based  on  the  work 
of  Havelock.2  The  first  classical  solutions  for  the  circular  cylinder 
(Lamb)17  and  the  sphere  (Havelock)18  have  contributed  much  to  the  general 
understanding  of  the  subject,  but  these  solutions  must  be  applied  with  great 
caution  to  problems  connected  with  elongated  bodies.  The  reason  herein  is 
the  extreme  sinplicity  of  the  cylinder  and  sphere;  the  resistance  curves  of 
these  bodies  do  not  show  the  characteristic  interference  effects  which  are 
peculiar  for  prolate  bodies  of  revolution.  From  physical  reasoning  we  infer 
at  once  that  in  the  latter  case  two  similarity  parameters  are  involved:  the 
common  Froude  number  F  =  V/VgL  referred  to  the  length  L  and  a  parameter 
characterising  the  depth  of  immersion  f,  say  f/L  or  the  depth  Froude  number 
Ff  =  U/VgT,  while  the  shape  of  the  wave- re si stance  curves  for  the  circular 
cylinder  and  the  sphere  depend  only  upon  F^.,  and  the  parameter  f/L  appears 
as  a  scaling  factor  only.  Thus,  for  instance,  the  peak  of  the  resistance 
curve  is  located  at  Ff  =  1  for  the  cylinder  and  Just  below  Ff  =  1  for  the 
sphere.  It  can  be  easily  shown  that  this  unity  value  of  the  depth  Froude 
number  has  no  special  significance  for  the  wave  resistance  of  a  very  elon¬ 
gated  body  of  revolution. 

Solutions  for  the  spheroid  and  general  ellipsoid  due  to  Havelock3,4 
lead  to  results  which  admit  of  qualitative  and  even  of  quantitative  esti¬ 
mates  of  the  resistance  of  "normal"  bodies  of  revolution.  The  importance 
of  the  spheroid  for  general  research  on  the  subject  cannot  be  overemphasized. 

Using  Havelock’s  general  expression  valid  for  a  plane  source-sink 
distribution,3  formulas  were  obtained  which  represent  the  wave  resistance 
of  a  rather  wide  class  of  bodies  of  revolution.5  By  these  formulas  the 
resistance  of  various  forms  has  been  investigated;6  especially,  some  endea¬ 
vors  were  made  to  find  forms  of  least  wave  resistance.5  These  forms  vary 
obviously  with  the  Froude  number  and  to  a  lesser  degree  with  the  depth 
parameter  f/L.  The  rather  striking  results  found  in  this  way  were  checked 
experimentally  and  good  agreement  between  theory  and  measurements  was  es¬ 
tablished  as  to  the  general  trend.6 
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As  with  surface  vessels,  theoretical  forms  of  least  wave  resis¬ 
tance  are  symmetrical  with  respect  to  the  midship  section.  Any  departure 
from  symmetry  causes  an  increase  in  wave  resistance,  and  this  increase  can 
become  appreciable  in  some  ranges  of  Froude  numbers  when  the  asymmetry  is 
pronounced.  The  degree  of  asymmetry  can  be  described  in  the  usual  way, 
though  roughly,  by  the  location  of  the  center  of  buoyancy  xQ,  or  the  dif¬ 
ference  of  the  prismatic  coefficients  ^p,  ^  of  the  fore  and  afterbody. 

For  instance,  a  difference  ^p  -  s  0.2  means  a  large  deviation  from  sym¬ 
metry.  Again,  the  resistance  results  are  qualitatively  supported  by  experi¬ 
ments.6 

An  extensive  hydrodynamic  study  of  bodies  of  revolution  is  under¬ 
way  at  the  Taylor  Model  Basin.  It  is  based  on  a  systematic  variation  of 
analytically  defined  forms.2, X1  As  an  extension  of  this  work  it  was  decided 
to  make  a  more  comprehensive  theoretical  investigation  on  the  wave  resistance 
of  bodies  of  revolution.  This  is  the  subject  of  the  present  report. 

In  Section  1  of  this  report  polynomials  are  discussed  which  are 
suitable  for  the  representation  of  hydrodynamic  singularity  distributions 
(doublets,  sources  and  sinks);  to  the  first  approximation  the  equation  of 
the  doublet  distribution  coincides  with  the  equation  of  the  sectional-area 
curve  except  for  a  scale  factor.7,8  A  class  of  curves  is  selected  which  in¬ 
cludes  the  TMB  Series2  generalized  by  one  additional  arbitrary  parameter. 

For  this  family  a  set  of  auxiliary  integrals  covering  a  large  range  of 
Froude  numbers  has  been  tabulated.  The  values  of  these  integrals  furnish 
immediately  the  variable  part  of  the  wave  resistance  of  the  simplest  forms 
(parabolas  of  the  type  1  -  $n).  In  the  general  case  the  wave  resistance  is 
given  by  a  quadratic  form  of  the  parameters  of  the  body  in  which  the  tabu¬ 
lated  values  appear  as  coefficients.  Thus  the  computation  of  the  wave  re¬ 
sistance  involves  only  some  multiplications  and  an  algebraic  addition. 

The  auxiliary  integrals  mentioned  have  been  computed  by  the  Bureau 
of  Standards.  A  short  description  of  the  work  involved,  contributed  by  Mr. 
Blum  of  that  Bureau,  and  tables  of  functions  are  found  in  Appendices  II  and 
III. 

As  mentioned  before,  the  resistance  formula  for  a  line  distribution 
of  singularities  used  throughout  this  report  follows  immediately  from  a  more 
general  expression  due  to  Havelock3,5  and  therefore  will  be  called  Havelock's 
integral . 

Using  the  tables  annexed,  resistance  curves  are  plotted  for  vari¬ 
ous  basic  forms  of  sectional-area  curves  (doublet  distributions);  they  cover 
three  depths  of  immersion  ratios  f A  except  for  the  spheroid  where  a  fourth 
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t/L  ratio  has  been  added.  Special  Investigations  are  made  on  the  Influence 
of  asymmetry,  and  some  examples  of  resistance  curves  refer  to  forms  selected 
from  the  TMB  Series. 

Following  an  earlier  attempt  distributions  of  least  wave  resistance 
are  investigated.5  Former  results5  are  checked  and  refined.  Particularly, 
the  distributions  obtained  lead  to  rather  peculiar  "swan-neck"  forms,  for 
higher  Froude  numbers.  Finally  it  is  shown  how  systematic  sets  of  resis¬ 
tance  curves  can  be  obtained  for  families  of  sectional-area  curves  (doublet 
distributions ) . 


2.  THE  REPRESENTATION  OF  SINGULARITY  DISTRIBUTIONS 
AND  SECTIONAL-AREA  CURVES  BY  POLYNOMIALS 

2.1 .  CONNECTION  BETWEEN  BODY  FORM  AND  GENERATING  HYDRODYNAMIC  SINGULARITIES 

In  establishing  a  relationship  between  body  form  and  generating 
hydrodynamic  singularities  two  well-known  problems  can  be  formulated: 

a.  Given  a  distribution,  find  the  shape  of  the  body  (sectional-area 
curve  A(x)). 

b.  Oiven  a  body  form  (sectional-area  curve  A(x)),  establish  the  appro¬ 
priate  distribution. 

In  the  present  report  we  disregard  the  complications  connected 
with  problem  b  and  treat  it  in  a  very  approximate  way.  The  contemporary 
rudimentary  state  of  knowledge  on  problems  of  wave  resistance  Justifies  this 
procedure  to  some  extent;  our  investigation  deals  essentially  with  resistance 
properties  of  hydrodynamic  distributions  and  merely  some  assumptions  are  made 
as  to  the  probable  shape  of  the  bodies  generated  by  these  distributions. 

Thus  two  essential  sources  of  error  are  involved  when  investigating 
the  wave  resistance  of  bodies  of  revolution: 

a.  The  approximate  character  of  the  wave-resistance  theory,  and 

b.  The  generally  admitted  approximation  that  for  a  given  body  the 
deep -immersion  distribution  of  singularities  can  be  used  instead  of  the 
actual  distribution  valid  for  near-surrace  conditions. 

The  second  assumption  (b)  appears  to  be  a  serious  one  when  the 
body  is  close  to  the  surface.  It  has  been  proved  by  Havelock*  that  it  leads 
to  inconsistent  results  with  respect  to  added  masses;  however,  by  following 
numerous  comparisons  between  theoretical  and  experimental  results  referring 
to  surface  ships  it  works  reasonably  well  when  applied  to  the  resistance 
problem. 

In  thv?  present  report  the  assuiptlon  will  be  made  that  the  shape 
of  the  body  generated  by  singularities  moving  close  to  the  surface  le 
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identical  with  the  shape  of  the  corresponding  body  generated  by  the  same 
singularities  in  an  unbounded  fluid. 

It  is  well  known  that  in  the  latter  case  one  can  construct  the 
contour  of  a  body  of  revolution  for  any  given  singularity  distribution  along 
the  axis;  auxiliary  tables  for  this  work  are  available,2, 10  especially  for 
cases  in  which  the  distribution  is  given  by  polynomials.  Plat  noses -as 
discussed  by  Weinstein14  -  will  not  be  dealt  with  in  the  present  report, 
although  it  is  possible  that  such  forms  are  advantageous  from  a  point  of 
view  of  wave  resistance  at  high  Proude  numbers.  When  dealing  with  "normal" 
shapes,  the  important  approximation  developed  by  Weinig7  and  Munk*  holds; 
i.e. ,  for  very  elongated  bodies  the  sectional-area  curve  of  the  generating 
body  A(x)  is  affine  to  the  doublet  distribution  /i(x).  This  approximation 
will  be  used  throughout  the  present  report  although  its  limitations  should 
not  be  forgotten. 

Some  explanation- if  not  definition -must  be  given  as  to  the  concept 
of  a  "normal"  shape  of  a  doublet-distribution  or  a  sectional-area  curve.  It 
means  essentially  a  curve  whose  trend  is  similar  to  sectional-area  curves  of 
common  ocean-going  ships;  these  curves  generally  are  monotonic  with  not  more 
than  one  point  of  inflection  in  the  fore  and  afterbody. 

Since  for  closed  bodies  the  source-sink  distribution  c  (x)  is  the 
derivative  of  the  doublet  distribution  u(x)  the  latter  is  monotonic  over  the 
range  of  the  forebody  when  a(x)  consists  only  of  sources  in  the  same  range. 
This  condition  (though  not  necessarily  a  required  one)  is  sufficient  to  ob¬ 
tain  bodies  such  that  the  circle  of  curvature  at  the  nose  lies  inside  of  the 
meridian  contour. 

We  mention  some  conditions  under  which  the  affinity  between  the 
doublet  and  the  sectional-area  curve  becomes  strained: 

a.  Por  larger  values  of  the  elongation  D/L  the  divergence  between  the 
sectional-* rea  curve  A (x )  and  the  doublet  distribution  u(x)  becomes  more 
pronounced  even  for  '‘normal"  shapes.  This  divergence  can  be  roughly  de¬ 
scribed.  First ,  in  the  mutual  relation  of  the  prismatic  (area)  coefficients 
which  are  the  decisive  form  parameters  of  the  two  curves- the  one,4d>  de¬ 
noting  the  prismatic  or  area  coefficient  of  the  distribution,  and  the  other, 

,  the  corresponding  one  for  the  sectional -area  curve-the  following  state- 
s 

ment  holds  for  a  wide  class  of  normal  bodies:9, 19 
for  finite  D/L 

>  *d  when  *d  <  2/3 
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*B  <  *d  when  *d  >  2/3 
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Figure  2  -  Spheroid.  Sectional-Area 
Curve  A,  Doublet  and  Source-Sink 
Distribution 


The  equality  ^  Is  valid  only 
for  the  ellipsoid;  see  Figure  2. 
Second.  In  the  prismatics  a  differ¬ 
ence  arises  between  the  length  of 
the  body  L  and  the  distribution  2a, 
2a  being  smaller  than  L.  For  the 
spheroid  the  relative  difference 

tu  L-2a  ~  D2  s  b2 
_55_=2L2  2a2 


where  f  depends  on  the  shape  of  the  distribution,  especially  at  the  ends. 
(Since  this  problem  is  being  thoroughly  Investigated  by  L.  Landweber  of  the 
Taylor  Model  Basin,  we  confine  ourselves  to  these  brief  remarks.) 

b.  When  conplicated  "abnormal"  distributions  like  "swan  necks"  or 
curves  with  very  steep  ends  are  investigated  (for  Instance,  Ranklne's  ovoid) 
the  divergence  between  these  distributions  and  the  sectional -area  curve  can 
become  appreciable  even  for  smaller  D/t. 


2.2.  REPRESENTATION  BY  POLYNOMIALS 
2.2.1.  General  Remarks 

In  former  reports  polynomials  have  beer  used  for  the  representation 
of  the  generating  doublet  (source  and  sink)  distribution  along  the  axis.5,,,1#,u 

The  doublet  and  source-sink  distributions  n  (x) ,  <r(x)  can  be  split 
up  Into  dimensional  factors  nQt  oQ  and  variable  dimensionless  parts  /i*U), 

**(():  »(x)  ■  »</*(!) 

o(x)  =  oQo*(t) 

with  $  =  x/a;  see  Figure  3b. 

The  dimensional  factors  will  be  established  later;  In  the  succeed¬ 
ing  discussion  the  functions  //*($)  and  <r#($)  will  be  treated  in  the  same 
way  as  ship  lines  and  their  derivatives.  Generally  following  Munk  and  Weinig 
the  doublet  distribution  #'*(()  Is  identified  with  the  sectional -area  curve 
A*U)  and  the  symbol  rj  is  used  for  both  of  them.  Actually  the  resistance 
computations  refer  to  given  distributions  for  which  the  corresponding 
sectional-area  curves  can  be  easily  calculated2’ lc  when  Munk's  approximation 
is  not  accurate  enough- as  for  Instance  in  cases  dealt  with  in  Section  5. 

The  first  adequate  representation  of  ship  lines  by  polynomials  is 
due  to  Taylor;11*12  the  equations  obtained  are,  however,  suitable  for  a 
separate  description  of  the  fore  or  afterbody  only.  Taylor  locates  the 
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Figure  3a  -  Landweber's  Axes  Figure  3b  -  Present  Axes 

Figure  3  -  Systems  of  Axes 


origin  at  the  bow  or  stem.  The  present  writer  has  proposed13, 14  other  sets 
of  polynomials  referred  to  a  system  of  axes  with  an  origin  located  midships. 
This  approach  has  definite  advantages  when  investigating  the  wave  resistance. 

Landweber2  has  generalized  Taylor's  equation  by  adding  one  more 
term  and  by  introducing  appropriate  boundary  conditions;  he  uses  the  ex¬ 
pression  obtained  as  the  equation  of  the  sectional-area  curve  of  a  four- 
parameter  form.2  The  parameters  are  interpreted  geometrically  as  the  pris¬ 
matic  coefficient,  the  location  of  the  maximum  section  along  the  axis  and 
the  nose  and  tall  radii  of  curvature.  It  will  be  immediately  shown  that 
Landweber's  equation  transferred  to  an  origin  at  the  midship  section  can  be 
split  up  into  a  two-parameter  symmetrical  and  a  two-parameter  skew  part  with 
respect  to  this  section;  thus  expressions  are  obtained  for  which  the  wave 
resistance  can  be  calculated  in  a  simple  way. 

2.2.2.  The  TMB  (Landweber)  Class  of  Bodies  and  Some  Generalizations 

The  TMB  (Landweber)  class  of  bodies  of  revolution  is  given  by  the 
equation  of  the  sectional-area  curve 

y2  «  a'x  +  a'x2  +  a'x3  +  a'x4  +  a'x5  +  a'x6  [1] 

11  21  31  41  51  61 

referred  to  axes,  as  shown  in  Figure  3-  We  transform  the  equation  of  the 
body  by  shifting  the  origin  to  the  midship  section  x  =  0.5»  reversing  the 
direction  of  the  axes,  and  putting  the  length  of  the  body  equal  to  2. 

Thus  for 

x  =  0  £  =  +  1 

i  ' 

xx  =  0.5  *  =  0 

x  *=  +  1  t  «  -  1 

l 
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The  transformation  is  given  by 


-2(xj -  0.5)  or 


The  resulting  equation  Is 


y2  ■  A  +  A  (  +  A  **  +  A  t*  +  A  +  A  *5  +  A  ** 
0  1'  2'  3'  4  s  •' 


where 


A  -  -/  nlnjLLlln.r_ 

3  <*s  ■>  .  <41 


n(n  -  1)  *n 
-Ll  etc. 


3  1  2x3x2" 

Equation  [3]  can  be  spilt  up  Into  a  symmetrical  and  an  ant  lsymmetr  leal  part 


y2  -  Y  -  i„  +  t# 


Y.  -  A„  +  A  {*  +  A  (*  +  A.<* 

Yn  «  A  i  +  A  4®  +  A  e5 
a  1'  3'  5' 


The  obtained  fora  [3]  has  definite  advantages  when  calculating  the  wave  re¬ 
sistance  since  the  latter  Is  the  sun  of  the  wave  resistance  corresponding  to 
the  synetrlcal  and  antlsymetrlcal  part  computed  Independently. 

Going  further,  we  derive  from  [3 1  the  following  simple  properties 
of  the  Landweber  bodies: 


vn 


The  coefficient  kQ  can  be  factored  out  and  merged  into  a  dimensional  constant 
which  defines  the  midship  section.  Thus,  the  normal  form  of  our  polynomial 
Is  obtained 


”  =  1  -2\<" 


with  a<  -  — 

1  Ao 

The  symmetrical  part  of  [4c]  Is  a  two-parameter  family 
>7flU)  -  1  -  a  j2  -  a J*  -  a  =  1  -  e#  -  •  («a-  *‘) 

B  2  4  6  2 

because  from  the  boundary  condition 

»,(' )  ■  0 

a  *>  1  .  a  -  a 

<  2  4 
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Such  families  have  been  called  "basic  forms"  by  the  present  writer1*  and 
designated  by  (2,4,6;$;t)  since  the  arbitrary  parameters  a  a  can  be  de- 

rl  2  4 

termined  by  the  prismatic  coefficient  4  “I  and  by  Taylor's  tangent 

value  t  *  -  dijO  )/$$. 

It  is  thought  that  the  Landweber  Series  [1 ]  meets  almost  all 
reasonable  requirements  as  to  wave -resistance  properties  presented  by  prac¬ 
tice  although  only  two  arbitrary  parameters  ^,t  are  at  our  disposal  for  the 
main  symmetric  part.  The  reason  for  this  assumption  is  that  from  investi¬ 
gations  on  surface  ships  it  is  well  known  that  area  curves  of  fine  ships, 
based  on  the  basic  family  equation  (2,4,6;^;t)  are  advantageous  in  the  range 
of  high  and  medium  Proude  numbers.  At  low  Froude  numbers  other  polynomials 
are  preferable  but  there  the  wave  resistance  of  submerged  bodies  becomes 
rather  negligible. 

We  have,  however,  introduced  an  additional  term  ag£*  for  which 
auxiliary  wave-resistance  functions  are  also  tabulated  in  this  report;  thus 
more  elaborate  investigations  can  be  performed  using  the  polynomial 


”8  ■  ’  'I 


1 


The  asymmetric  (skew)  part  is  the  function 


n  =  a  (  +  a  {*  +  a  {s 
a  is  s' 

factoring  out  we  write 


[4d] 


tjq  =  a  it  ~  a  U  +  b  +  b  $5)  [4e] 

a  i  'a  i  ’  3  5 

Obviously  the  resultant  curve  17  =  rj_  +  rja  can  have  its  maximum  section  out- 

5  a, 

side  of  (  =  0  and  the  area  of  this  section  will  generally  differ  from  one. 
This  slight  complication  does  not  involve  any  difficulties  in  actual  work. 
Let  us  investigate 

*?;  =  I  +  b3*3  +  b5*5  [4f  3 

This  trinomial  has  to  comply  with  the  conditions 

hj(o)  =  0 

nj(+!)  =  ij(-n  =  0 

whence 

b  =  -(1  +  b  ) 

5  3 


t 

I 

\ 


i 

i 

$ 
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thus 

rtl  =  *  +  b3*3  -  0+b^*5  [4g] 


The  only  arbitrary  parameter  b  can  be  fixed  by  one  additional  condition;  as 
such  we  choose  the  tangent  value  t*  at  the  bow  (at  the  stem  the  correspond- 

a 

ing  value  is  -t*) 

a 


3i7jO) 

d* 


=  +  4  +  2b 

3 


hence 

b  =  -  2  +  t*/2  [4h] 

3  a 

the  corresponding  tangent  value  ta  of  17  ,  Equation  [4e],  is  obviously 

a  a 

t  =  a  t* 
a  x  a 

The  table  below  shows  some  examples  of  skew  forms.  The  parameter 
4*  =j0  4  *  Is  an  area  coefficient  referred  to  the  unit  square.  Plots  of 

$  -  $3  ,  $  -  £ 5  and  some  other  "skew"  forms  used  in  the  TMB  Series  are  shown  on 
Figure  5-  The  actual  skew  part  17  contains  additionally  the  "strength  para- 

d 

meter"  a  ;  see  Equation  [4e]. 


0  1 

2 

4 

"a* 

*0  -  t2)2  t  -  1-5*3  +  O.5*5 

*  -  *3 

«  -*5 

drj£/d{ 

1  -  6$2  +  5$4  l  -  4.5$2  +  2.5*4 

1  -  3*2 

i  -  5*4 

<A* 

1/6  =  0.166...  5/24  =  0.2083... 

1/4 

1/3 

Our  numeric  evaluations  are  primarily  based  on  Equations  [4c]  and 
[5] -which  are  stated  below-but  the  theoretical  treatment  will  be  carried 
out  along  more  general  lines. 

Extended  investigations  have  been  made  by  Landweber  and  Gertler2 
on  the  influence  of  an  additional  term  ajx7  on  the  form  of  the  body  when  the 
geometric  parameters  are  kept  constant. 

Using  our  system  of  axes  it  is  easy  to  perform  similar  investiga¬ 
tions  for  the  symmetric  and  asymmetric  part  of  [4c] 


1 


Figure  4  -  Dimensionless  Sectional-Area  Curves  A*(£ ) (Doublet  Distribution^ 
M*\0)  of  Some  Simple  Bodies  Symmetric  with  Respect  to  the  Midship  Section 


Figure  5  -  Examples  of  the  Antisymmetrical  (Skew)  Parts  of  Sectional -Area 
Curves  A*U)  (Doublet  Distributions  */*(£))  Belonging 

to  the  Family  £  +  b  -  (1  +  b  ){5 


12 


By  adding  to  ijq  terns  with  arbitrary  parameters  a7  and  a,  ,  a  manifold 
ijx  =  +  aT$7  +  a$$*  Is  obtained. 

The  polynonlal  [4c]  Is  completely  defined  by  the  four  geometrical 

parameters 


U  # 


a»  (l ) 


pressed  by 


where 


2)  ts  s  '  d< 
dr,JO) 

“ft—*! 

au.d) 

4)  -fr-s -  'V?  -  -fca 


When  ij  has  to  comply  with  the  four  equations  [5]  It  can  be  ex- 


»ix  “  »j0  +  cs43»iU)  +  cad4»?(0 

*,»<«>  =  -  5«*  +  7t '  -  3« 


complies  with  the  conditions 


=  0; 


satisfies 


,  diMO)  8Am n(l  ) 

AMO)  =  4,11(1)  - i - 0 


=  *3o  -  t2y 


dA  n(0)  ad  n(l) 

V(0)  '  V(,)  =  ^  *  0 


Thus,  an  addition  of  the  functions  A  ,  A  to  n  does  not  influence  the 

34  0 


boundary  conditions,  [5].  The  shape  of  the  curves  A^rjl^)  and  A  »?U)  Is 
shown  In  Figure  13-  The  advantage  of  this  representation  Is  obvious. 
While  In  the  equations 


’  -  "8  +  "a 


»*U)  -»*(«)  +«*(«) 
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the  symmetrical  (even)  terms  17  ,  ft?  are  the  main  parts,  obviously  in 

5  s 

an  _  d's  ,  **a 
3*  '  3*  a* 
or 

*•(«)  *  *JU)  +  **(«) 

the  odd  terms  d>j_/3|,  <r*  become  the  main  part. 

S  S 

2.3  CONNECTION  BETWEEN  STRENGTH  OP  SINGULARITIES  AND  BODY  SHAPE 

The  next  consideration  is  to  establish  the  dimension  factors  ft  ^ 

o 

and  oQ.  The  flux  through  the  midship  section  may  be  written  as 

Q  -  C(|,  /i*)*^  [7] 

Here  the  coefficient  C(b/a,  /1*)  is,  as  indicated,  a  function  of  the  elonga¬ 
tion  ratio  b/a  s  d/L  and  of  the  shape  of  the  distributional*.  For  very’  large 
elongations  C(b/a,  /i*)-^1,  but  for  shapes  and  values  b/a  used  in  actual  op¬ 
eration  C  differs  from  one. 

A  closer  investigation  of  the  coefficient  C  will  be  given  elsewhere 
by  L.  Landweber;  for  the  present  purpose  we  Introduce  C  as  a  correction  fac¬ 
tor  which  improves  the  accuracy  of  Munk's  or  Weinig's  approximate  affinity 
theorem  mentioned  on  page  5.  The  dependence  of  C  upon  ft*,  although  apparent¬ 
ly  negligible  within  the  range  of  presently  used  submarine  hull  forms,  shows 
some  interesting  features.  Earlier  brief  Investigations  lead  to  the  follow¬ 
ing  table  for  C(b/a,  1 1 *)  (Reference  5). 


b/a  *  D/li 

1/4 

1/6 

1/8 

1/10 

0 

#»*(«) 

*d 

C(b/a,#i*) 

(1  -<2)2 

0.533 

1 .172 

1.093 

1 .060 

1.043 

1 

1  -*2 

0.660 

1 .192 

1 .092 

1 .054 

1 .036 

1 

1  -  3.0825**+  0.165*10  +  1.9175*12 

0.820 

- 

1 .01 24 

1.008 

1 

Prom  these  results  we  gather  that  C(b/a,  ft*)  values  for  normal 
submarine  shapes  apparently  can  be  estimated  from  the  spheroid;  an  empirical 
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formula  C( b/a)  s  i  +  3b2/a2  may  fit  the  facts  reasonably  well.  For  fuller 
bodies  lower  values  seem  to  be  suitable. 


The  constant 


is  therefore  obtained  as 


"o  ■ 


*o  =  lc< }’  !7a] 

The  flux  [73  or  the  strength  of  the  doublet  distribution  at  the 
midship  section  must  be  somewhat  higher  than  the  product  of  the  cross  section 
times  the  speed  of  advance. 

For  the  source ,  we  have 

•o  ■  f  c<r  "*)  rB  18) 

3.  EVALUATION  OF  HAVELOCK'S  INTEGRAL 
3.1 GENERAL  CONSIDERATIONS 

The  save  resistance  experienced  by  a  continuous  doublet  sheet  p, 
distributed  over  a  vertical  plane  and  moving  uniformly  on  a  straight  hori¬ 
zontal  path,  has  been  calculated  by  Havelock.3  Concentrating  the  distribu¬ 
tion  fi(x)  along  a  horizontal  straight  line  we  obtain  immediately 

R  -  i6ttPkJ  f**  |P2  +  Q2J  sec50d0;  KQ  =  [9] 


with 


/  +  a 

p(x)  cos  (Kqx  sec  0)dx  =  exp(-KQf  sec2#^  [9a] 

r+a 

Qx  =  exp(-KQf  sec z9) J  /i(x)  sin  (KQx  sec  0)dx  =  exp(-KQf  sec2#^  [9b] 


hence 


R  =  l6»rpK£  j  (p2  +  q2)  exp(-2KQf  sec20)  sec50d0 


Using  a  source -sink  distribution  we  obtain  similarly 


R  =  1 6np  K2  f  /2  (p2  +  q2)  exp(-2K0f  sec20)  sec30d«  [10] 
J  0 

r  +  a 

p  =  a[x)  sin  (K  x  sec  0)dx  [10a] 

J  -a 

[*• 

q  ■  o(x)  cos  (K.x  sec  9)d9  [10b] 

J  -  A  O 


[10b] 
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Introducing  dimensionless  coordinates  x  -  a{  and  the  expressions 

m(x)  - 

o(x)  •  <rQo*(t) 

various  forms  of  the  integral  for  R  can  be  derived  for  purposes  of  numerical 
evaluation. 

We  confine  ourselves  to  the  source-sink  integral.10  Splitting  up 
into  a  main  antisymmetrical  and  a  symnetrical  part 


*#U)  -<rJU)  +  **(0 


and  remembering  that  an  integral  taken  over  an  odd  integrand  between  limits 
of  equal  absolute  value  but  opposite  sign  vanishes,  we  obtain  with  the  desig¬ 
nation 


*o 


V 


U2  2F2 


R  =  4C2wpg  j  /2  exp [-4  £  yQ  see2#]  sec3fl.  [11  ] 

{ //?'«  sin<ro  <sec  #)d«}2  008  <V sec  «)d*J2jd*  * 

const  J  2  exp  [-4  £  yQ  sec2flj  sec30[p#2  +  q*2]dfl  [12] 
We  Introduce  further  polynomials  for 

u*U)  -  1  -£****  03] 

hence 

•*(«)  -  ^Tn  anr“  [14] 

or 

**U)  =  *JU)  +  **(«)  »  -2’2ka2Jk<  2k_1  (2»H  (l4al 

k  * 

with  k,  m  as  integers. 

For  the  main  antisymmetrical  part  the  intermediate  integral  p#  becomes 

P*  *  J  «J(?)  Sin  (r0?  sec«)d(  *  -^2k  a(  J"  i**'lsln  (r0<  sec*)d{ 

=  -T?k  H2l_i(y0  sec  8)  [15] 

with 

M2»-i*>’o  sec  «)  -  Mjt  i(r)  -  J1  <”  "l  Bln(r0t  sec  a)d{  «  J1{2*"lsin()’<)d{ 

[16] 


*  V 

* 


1 6 


Here  for  brevity  the  designation  y  =  yQ  sec  0  has  been  introduced,  [ 1 6 ] . 
For  the  symmetrical  (even)  part 

q*  c  f  <7?U)  cos(y t  sec  0  )d£  =  -  y"(2nH-l  )a  M 1  ( y  sec  0  )  (17) 

Jc  S  O  ■  m  <  2m  +  1  2m  O 

with 

M'  (yrt  sec  0 )  =  M'  (y)  -  M1  =  f  f2m  cos(y  (  sec  0 )d{  =  f  /  2m  cos  y  d{ 

2m  O  2m  2m  !  O  J 

Jo  Jo  (l8] 

inserting  [15]  and  [l6j  into  [12]  one  obtains 

R  -  4C2  npg  l'  ,0  J’/’exp  [-4  £  y0  sec2  «]  ((^>  2  + 

+  ( ^7  (2mt1  )a  M'  )2l  sec30d0  [19] 

This  formula  is  suitable  for  numerical  computations  above  in  special  cases, 
since  tables  of  the  functions  M  (y),  M*  (y)  are  available  and  will  be 

24“ 1  2m 

published  in  a  TMB  Report. 


3.2.  TABULATION  OF  RESISTANCE  INTEGRALS  FOR  A 
FIVE-PARAMETER  CLASS  OF  BODIES 

As  mentioned  before,  auxiliary  integrals  have  been  prepared  for  the 
three-parameter  symmetric  distributions  of  Equation  (5],  u*(()  (asymmetric 
inaJU)). 


*:«>  - 1  -Z  -n 


a  and  a 


s' T '  ~n  ’  -  '  a^  “  "  mn  * 

2,4,6,* 

and  the  one-parameter  skew  distribution  //*({)  (symmetric  in  a|U)) 


M*U)  -  *  +  b<8  -  (1+b  )*5  and  a*U)  »  1  +  3b  *2  -  50+b  X4 
a  s  3  a  33 

The  computations  are  based  on  a  slightly  different  form  of  R  (see  Appendix 
II).  Substituting 


y  =  yQ  sec  0 


sec  0  =  y/y  tg  0  =  ^Ty/y  )2  -1 


one  obtains 


dy 


yQ  sec2  0  sin  0  d0; 


yQ  seca  0d0  = 


(y/y0)2dy 


[/(y7>'0 ) ; 


-1 


hence 


R  -  4C2  exp[-4  £  £]  {(Z^^fr))2  ♦ 

+  If <a,H  t»i 
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putting  for  abbreviation 


4C2  nfig  ~  =  const  =  C, 

(r/r0)2 

~r==  “  f(y) 

^(y/r0)2  -i 


R  -  C0  J*  exp[-  £  f(r)  {£21  2]  aal  «2.  H,(  + 

+  £(2r+1H2^1)a,„ialitlR'r  M'jdy  [21] 
with  1,  J,  s,  r  Integers.  R  can  be  built  up  of  terms  of  the  type 

f*exp  [- ^  £lf (r)H  „(r)R,  „(rWr  -  »  ...  ....  [22] 

Jyo  L  L  yQ J  2‘  1  2;  1  2‘  1,  2j  1 

for  the  symmetrical  part  of  the  sectional-area  curve 
and 

0  U 

for  the  skew  part.  The  final  result  Is  therefore  obtained  as  a  quadratic 
form  in  the  parameters  a  or,  better,  nan 

R  =  Crt{  T  21  2j  a  am  +  >T(2rf1 )  (2s+1  )a  a  TTp  l[24] 

OVf*  2.  2j  2 1  "I,  2)  "I  Tjf  *r+l  t*  +  l  2r,  2*  J 

with  the  tabulated  integral  values  T7?  as  main  parts  of  the 

2i  "1,  2j  -1  2r ,  2s 

coefficients  21  2J  77}  ,etc. 

2.  -  1,2>-X 

We  mention  again  the  fortunate  circumstance  that  the  contributions 
to  the  wave  resistance  due  to  the  symmetrical  and  antlsymnetrlcal  parts  can 
be  calculated  Independently  and  added. 

Returning  now  to  a  family  of  distribution  curves  given  by  Equation 
[4c]  but  generalized  by  one  additional  term  a8  f, : 

i  - 1  -  i\<"+  ■>,  t'-§f  -  -  £»  vn" +  k/  [41) 

The  wave  resistance  can  be  calculated  by  the  functions 


^11 

**10 

n, 

K 

<2  K4 

”>33 

^33 

™ 22 

^24 

^33 

^3, 

^44 

i9 


tabulated  In  the  Appendix  III.  The  integral  R  and  the  functions  7f\  and  7Jf 
depend  upon  the  two  parameters  yQ  ~  1/2P2  and  f/L.  The  tables  have  been  pre¬ 
pared  for  a  range  0.5  ^  Y0  £  10  and  f/L  =  0.125,  0.25,  0.50.  Additionally, 
for  >}  an  intermediate  depth  of  immersion  ratio  f/L  =  0.1875  has  been  intro¬ 
duced.  From  the  wave  resistance  integral  it  follows  immediately  that  the 
ratio  depth  of  immersion  over  length  f/L  is  theoretically  preferable  to  the 
more  commonly  used  ratio  f/D,  since  f/L  appears  explicitly  as  factor  of  the 
exponent  of  the  e-function  under  the  integral.  With  elongated  bodies  the 
ratio  b/a  or  D/L  influences  primarily  the  constant  CQ  =  4C2  npg  b4/a  only, 
though  in  a  very  decisive  way.  Although  the  lower  speed  limit  =  10 
(F  =  0.224)— up  to  which  the  auxiliary  integrals  have  been  computed— is 
rather  high,  it  is  thought  that  for  normal  hulls  with  «*<2/3  moving  at  greater 
depths  than  D,  the  wave  resistance  becomes  uninportant  when  F  <~0.224.  The 
low-speed  range  may,  however,  be  interesting  in  connection  with  other  research 
problems . 

In  principle  the  wave- re si stance  equation,  [24],  solves  the  problem 
for  any  sets  of  an  within  the  family  following  [4i],  Actually  since  the 
relative  error  of  the  tabulated  functions  is  approximately  0.0001,  a  loss  of 
accuracy  may  occur— when  the  coefficients  an  reach  high  absolute  values  with 
alternating  signs.  It  is  not  probable  that  difficulties  of  this  kind  will  be 
important  in  connection  with  submarine  work;  besides,  they  can  be  overcome 
to  some  extent  by  plotting  suitable  sinpler  resistance  curves  and  by  inter¬ 
polating. 


4.  REPRESENTATION  OF  RESISTANCE  CURVES 
4.1  .  THE  DIMENSION  FACTOR  CQ  AND  DIMENSIONLESS  REPRESENTATIONS 

The  dimension  factor  in  Equation  [20],  CQ  c  4C2  npg  b4/a,  has  a 
rather  unusual  form,  but  it  will  be  widely  used  throughout  this  report 
because  of  its  theoretical  merits  and  the  comparative  ease  with  which  it  can 
be  connected  with  more  familiar  expressions.  We  rewrite,  in  terms  of  the 
displacement  A  , 


C  =  *7rb2  ?apg— -  C2  =  A  —  C2 
0  tfa2  tfa2 


R  =  r  .  |-*SL  or?  -  r 
Co  0  4  2b2C2  4  0  *a2 


[25] 
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Hence  we  can  immediately  derive  the  resistance  per  unit  displacement  for  a 
given  b/a  and  shape  when  rQ  Is  known. 

The  Introduction  of  the  displacement  A  in  [20]  Is  open  to  objection 
since  so  far  we  have  not  distinguished  between  the  length  of  the  body  and  the 
distribution.  We  repeat  the  definitions: 

2a  Is  the  length  of  the  distribution  along  the  axis 
L  Is  the  length  of  the  generated  body 
2b  =  D  Is  the  diameter  of  the  generated  body 
Obviously  for  the-  displacement  of  the  body  we  must  use  L  *  21.  Then 


C  =  Unpg  C2  b4/a  =  A  — -b2*  [25a] 

0  #a2  l 

Further,  the  ratio  b //  =  D/L  is  technically  more  Important  than  b/a;  hence 


or 


R__£_  H  a 

A  2C2  b2  l 


_R  _  2C2/D\2  l 
A~  <t>  IlJ  a 


Later  we  shall  use  another  coefficient 


r 

i 


ro  _  R  J_  /L\  _a 

-T-a  2C,\d)  i 


[25b] 


[25c] 


[26] 


One  should  not,  however,  overestimate  the  Influence  of  the  length  correction. 
For  the  spheroid 

Ti ?  s  (~  *  i  ¥)  o  *  *  m 

l.e.,  influences  the  C2  correction  by  less  than  10  percent.  Further,  even 
the  Introduction  of  the  more  important  C  factor  does  not  lead  to  an  exhaustive 
correction  since  we  know  that  not  only  the  midship  section  but  the  whole 
trend  of  the  curves  changes  with  increasing  b/a.  Thus  within  the  limited 
accuracy  of  the  present  wave-resistance  theory  we  generally  can  put  //a  ss  1. 

It  Is  of  course  important  to  use  all  approximations  in  a  consistent  and 
clearly  defined  way,  so  that  fair  comparisons  can  be  made. 

We  note  particularly,  that  for  the  spheroid 


1  -  VC2 


SiiWWiiffi 
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For  comparison  with  experiments  the  coefficients  cH  referred  to  the  wetted 
surface  S  Is  advantageous. 

We  write 


c  = 


4nC*b2/ b 


p/2  U2S 


o  sp2  la  f 


or  Introducing  a  surface  coefficient  S  (Reference  17) 


(271 


Cg  =  S/wDL 


■  _  Cf/b\*  b  1  ~ 

cw  “  ro  Cg  la/  C  pj  * 


[28] 


with  yQ  =  1/2P2.  For  elongated  spheroids  Cg  *  0.79* 

The  Importance  of  the  resistance  coefficient  cw  referred  to  the 
wetted  surface  S  Justifies  a  short  digression  on  the  calculation  of  S  for 
bodies  of  revolution.  Solutions  of  the  exact  expression  (Equation  [29])  can 
be  obtained  In  a  closed  form  In  exceptional  cases  only,  as  for  the  spheroid. 
Of  course  it  presents  no  difficulties  to  evaluate  the  Integral  numerically, 
but  a  simple  approximate  formula  can  be  derived  at  least  for  the  surface  area 
of  a  restricted  class  of  very  elongated  bodies  of  revolutions  complying  with 
the  condition  that  the  end  tangents  of  their  merldianal  contour  do  not  become 
vertical;  It  Is  similar  to  the  well-known  expression  for  the  length  of  a 
slightly  curved  arc,  see  Appendix  I. 


4.2  RESISTANCE  CURVES  OF  SIMPLE  SYMMETRICAL  BODIES 

Since  the  presentation  and  the  discussion  of  resistance  curves  is 
the  main  subject  of  the  present  report,  various  sets  cf  such  curves  have  been 
computed.  Essentially,  the  resistance  properties  of  the  following  three 
groups  of  body  forms  (distributions)  have  been  Investigated; 

(a)  A  set  embracing  a  wide  range  of  prismatic  coefficients,  which  fur¬ 
nishes  a  general  review  of  the  resistance  as  function  of  the  form  (IV, 2). 

(b)  A  set  dealing  with  four  TMB  models.  This  raises  the  problem  of 
the  Influence  of  asymmetry  with  respect  to  the  midship  section  ( IV, 3 ) . 

(c)  A  group  consisting  of  systematically  chosen  forms  belonging  to  the 
two-parameter  family  (2,  4,  6;  0;  t)  (VI);  for  the  same  family  some  calcula¬ 
tions  of  shapes  of  least  resistance  &re  presented  (V). 

The  procedure  adopted  leads  to  repetitions  which,  having  in  view 
the  Importance  of  the  subject,  have  been  thought  to  be  advisable.  Because 
of  the  complicated  dependencies  Involved  the  interested  reader  can  more 
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Figure  6  -  Wave-Resistance  Coefficients  r_  =  -  of  Symmetrical 

0  47rC2pg  bVa 

Bodies  as  Defined  on  Figure  4,  f/L  =  0.125 

easily  draw  conclusions  from  the  rather  comprehensive  plots  than  from  any 
text. 

We  are  mainly  interested  in  the  range  of  Froude  numbers  F  below 
and  at  the  maximum  of  the  large  hump  in  the  resistance  curve;  see,  for  in¬ 
stance,  Figures  6,  7  and  8.  Above  the  maximum  the  absolute  value  of  wave 
resistance  decreases  comparatively  slowly  with  growing  F,  but  the  ratio  wave 
resistance  to  frictional  resistance  drops  quickly.  Therefore,  at  high  speeds 


22 


Plgure  7  -  Wave- Resistance  Coefficients  of  Symmetrical 
Bodies  as  Defined  on  Figure  4,  f/L  =  0.25 

the  wave  resistance  of  elongated  bodies  such  as  torpedoes  represents  only  a 
small  part  of  the  total  drag.  It  has  been  shown  in  References  4  and  5  that 
in  the  limit  of  very  large  Froude  numbers  the  wave  resistance  becomes  pro¬ 
portional  to  the  square  of  the  displacement  or  rQ  to  ^2. 

In  general,  throughout  the  present  report  calculations  have  been 
extended  to  F  »  1  (yQ  =  0.5).  and  to  F  s  1.58  (yQ  =  0.2)  for  the  parabolic 
distributions  1  -  1  -  |4,  1  -  f*only.  From  an  approximate  Investigation 
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Figure  8  -  Wave-Resistance  Coefficients  of  Symmetrical 
Bodies  as  Defined  on  Figure  4,  f/L  =0.5 


it  appears  that  the  resistance  curves  R  (yQ)  plotted  over  yQ  have  a  vertical 
tangent  at  yQ  -  0,  but  no  attempt  has  been  made  to  draw  accurately  the  range 
of  curves  below  yQ  -  0.5. 

To  obtain  a  general  idea  of  the  wave  resistance  for  various  symme¬ 
tric  distributions  u(t)  (sectional-area  curves  A*($))  graphs  have  been  plot¬ 
ted  for  following  simple  cases:* 

*Aa  before,  by  syn»>etry  we  ween  synaetry  with  respect  to  the  mideection. 
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Figure  4  shows  these  sectional-ar#a  curves  and  Figures  6,  7  and 
8  the  corresponding  resistance  coefficients  as  functions  of  yQ  =  1/2F2,  with 
an  additional  non-equidlstant  scale  for  F.  The  choice  of  yQ  as  independent 
variable  yields  an  appropriate  picture  of  the  wave-resistance  values  at  high 
speeds . 

From  tb?  Figures  6,  7  and  8  a  rather  complete  understanding  of  the 
wave -resistance  properties  of  various  synnetrical  forms  can  be  derived!  Ref¬ 
erence  is  also  made  to  Figure  12  and  the  pertaining  discussions  in  the  text. 
The  influence  of  the  depths  of  immersion  follows  immediately  from  a  comparison 
of  Figures  6  through  8;  also,  cross  curves  can  be  plotted  over  f/L  as  the 
Independent  variable.  Figure  11  shows  this  dependency  for  4?^,  which  is 
the  resistance  function  of  a  spheroid  A*({)  «  1  -vf2,  with  yQ  =  1/2F2  as 
parameter.  We  note  that  with  increasing  depth  the  resistance  drops  more 
quickly  it  small  than  large  Proude  numbers  F.  This  is  rather  obvious; 
it  will  be  discussed  later  more  thoroughly  that  the  most  indicative  parameter 
is  the  ratio  f/x  ,  where  X  the  length  of  the  free  wave  is  X  =  2ttF2L. 

In  Figures  9  and  10  the  resistance  curves  for  three  depths  of  im¬ 
mersion  have  been  reduced  to  approximately  the  same  maximum  ordinates.  This 
rather  artificial  approach  yields  a  clear  idea  about  the  shift  of  the  last 
hump  (of  its  steep  rise  as  well  as  of  the  position  of  its  maximum)  to  higher 
Froude  numbers  with  increasing  depth  of  immersion;  it  further  emphasizes 
again  that  the  rate  of  decay  of  the  wave  resistance  with  increasing  depth  is 
much  higher  for  low  Froude  numbers  than  for  high  ones. 

Figure  12  represents  a  coefficient  r  *  R /A  a2/2C2b2  *  rQ/<*.  For 
approximately  constant  C2  (very  elongated  bodies)  and  given  a2/b2  ratio, 
r  ~  R /A,  i.e.,  the  figure  yields  a  comparison  of  the  resistance  per  unit 
displacement  for  various  forms. 

The  discussion  of  the  various  graphs  leads  to  the  following  summary 

results: 
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Figure  Q  -  Comparison  of  the  Shape  of  Wave-Resistance 
Curves  for  the  Spheroid  tj(4)  =  1  -  the  Curves 
are  Reduced  to  Approximately  Equal  Maxima 

A.  Small  depths  of  Immersion 

1 )  Within  reasonable  limits,  the  peak  value  of  the  R/d  curve  does  not 
depend  too  much  on  the  shape  of  the  body,*  especially  upon  the  prismatic  co¬ 
efficient. 

2)  The  merits  of  full  forms,  over  a  wide  and  possibly  Important  range 

of  Froude  numbers  0.35  ^  ^  0-50,  are  clearly  emphasized,  as  well  as 

3)  The  heavy  penalty  wnich  has  to  be  paid  for  high  prlsmatlcs  at 
lower  F. 


•If  oore  elaborate  results  are  desired  they  can  be  derived  from  Figures  28  through  35. 


?6 


Figure  10  -  Comparison  of  the  Shape  of  Wave-Resistance 
Curves  for  rj  *  (1  -  (2)z  Reduced  as  by  Figure  9 


B.  For  larger  depths  of  immersion  the  dependence  of  the  peak  values  of 
R/4  upon  <t>  becomes  more  pronounced;  the  advantage  of  high  prismatics  in  the 
range  mentioned  in  A ( 2 )  is,  on  the  average,  reduced. 

4.3.  RESISTANCE  CURVES  OF  ASYMMETRICAL  BODIES 

Further  curves  representing  the  wave-resistance  coefficients  of  the 
four  TMB  models  represented  in  Figures  14  and  1'5  are  shown  in  Figures  l6,  17 
and  l8.  Before  discussing  these  particular  asymmetric  models,  hov/ever,  an 
investigation  must  be  made  of  the  influence  of  asymmetry  on  the  resistance. 

Figure  5  represents  examples  of  asymmetrical  lines  belonging  to 
the  family  >??={+b{3-(l+b){5,  Equation  [4g], 
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Figure  11  -  Wave -Re si stance  Coefficients  of  the  Spheroid 
as  Functions  of  f/L  with  yQ  «  1/2F2  as  Parameter 
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Figure  12  -  Wave-Resistance  Coefficients  r  *  -£■  =  — - 

i  T  d  2Cab2 

of  Symmetrical  Bodies  Shown  on  Figure  4,  f/L  »  0.125 
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Figure  13  -  Distribution  Functions  Following  Equations  [6a]  and  [6b] 


The  curves  I  to  IV  have  been  derived  from  the  TMB  models  (Figure.' 

cl  cl 


14  and  15)  by  reducing  the  coefficient  of  (  to  unity.  The  procedure  of  ob¬ 
taining  the  symmetric  and  the  skew  part  from  graphs  is  obvious:  The  first 
one  is  the  arithmetic  mean  of  the  fore  and  afterbody  ordinates  rjm  - 
and  the  latter  one  the  difference  —  or  respectively. 

The  computation  of  the  wave  resistance  due  to  asymmetry  is  based 
on  Equation  [24]: 

For  the  trinomial 


a  (  t  +  b*3  +  b  *s) 

1  3  5 


8450 

8450 
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Figure  l6  -  Total  Wave -Resistance  Coefficients  and  Coefficients 
Due  to  Asymmetry  of  the  Four  TMB  Models  Shown 
in  Figure  14  and  Figure  1 5 »  f/L  =  0.125 
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Figure  17  -  Total  Wave-Resistance 
Coefficients  rQ  of  the  TMB 

Models  Shown  in  Figures 
14  and  15,  f/L  =  0.25 


Figure  18  -  Total  Wave -Resistance 
Coefficients  rQ  of  the  TMB 

Models  Shown  in  Figures 
14  and  15,  f/L  =  0.5 


with  the  derivative 


drj 

a 

d* 


-  ax(l  +  3y2  +  5b.*4) 


we  obtain 

R  =  C  a2  f77?'  +  9b2  71?'  +  25b2  W  +  6b  Iff  +  ICb  7 Tj*  +}0b  b  ttj'  1 

a  0  l  [  00  3  22  5  44  3  02  5  04  3  5  24J 

(33] 

Figure  19  shows  the  functions 

p 
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Figure  19  -  Wave-Resistance  Coefficients  rM  Due  to 
Antisymmetrical  Distributions  Following  Figure  5 


corresponding  to- the  distributions  -  $3,  £  -  £3  and  I*  shown  in  Figure  5, 

a 

where  curve  I*  is  derived  from  the  TMB  body,  Figure  14. 

The  "amount"  of  asymmetry  which  corresponds  to  the  equation  rj.  - 

a 

is  very  large,  but  by  assuming  the  strength  parameter  a^  <  1  (Equation  [4e]) 
more  usual  distributions  are  reached;  for  these  asymmetric  terms  the  wave- 
resistance  curves  are  obtained  simply  by  multiplying  the  ordinates  of  Figure 
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19  by  a*. 

The  resistance  curves  in  Figure  19  corresponding  to  f  -  f5  and  I| 
are  somewhat  similar  in  the  range  of  the  large  hump  and  the  ratios  of  their 
absolute  values  are  of  the  order  of  0-5.  In  the  range  of  the  second  hump  the 
ordinates  of  both  curves  are  small,  but  it  is  characteristic  that  here  a  much 
lower  resistance  corresponds  to  the  finer  line  I*,  rather  than  to  f  -  {5. 

o- 

We  return  now  to  the  four  TMB  models  designated  by  I,  II,  III,  IV 
shown  in  Figures  14  and  15-  In  these  figures  the  line  A*(f)  shows  the  sym- 
metrical  part  of  a  body.  The  resistance  results  are  plotted  on  Figures  l6,  17 
and  l8;#  in  them  the  lower  set  represents  the  contribution  due  to  antisymmetry 

y 

roa  “  knpgC*  b4/a  * 

the  upper  set  the  total  wave -resistance  coefficient 


The  computations  are  made  under  the  assumption  that  the  doublet  distribution 
A*(t).  With  the  model  number  rising  from  I  to  IV  the  prismatic  in¬ 
creases  and  the  asymmetry  decreases.  In  the  important  range  of  Froude  numbers 
0.50  >  F  >  0.35  the  finer  models  are  extremely  unfavorable  because  of  the 
low  prismatic  as  well  as  because  of  the  very  pronounced  asymmetry. 

When  comparing  the  total  resistance  values  a  slight  departure  from 
symmetry  generally  is  advantageous  because  of  viscous  effects.  It  has  also 
been  pointed  out  that  small  asymmetric  terms  do  not  increase  appreciably  the 
wave  resistance  even  in  the  most  sensitive  range  of  Froude  numbers,  say 
0.45  >  F  >  0.35;  this  is  well  supported  by  our  present  results,  for  instance 
by  Curve  IV.  Further,  the  obvious  fact  must  be  once  more  emphasized  that  an 
Immediate  comparison  between  symmetrical  and  asymmetrical  bodies— as  to  their 
wave-resistance  properties— is  only  feasible  when  the  sectional  area  of  the 
former  A*U)  is  the  even  part  of  the  sectional  area  of  the  latter 

A*U)  =  AJU)  +  AJU) 

It  is  entirely  possible  to  obtain  asymmetrical  forms  with  wave-resistance 
properties  which  are  superior  to  the  corresponding  ones  of  a  poorly  chosen 
symmetrical  form,  equal  prlsmatlcs  and  principal  dimensions  being  assumed. 

Similar  computations  have  been  performed  for  other  depths  of  im¬ 
mersion;  some  results  are  listed  in  Table  2  of  Appendix  III.  Obviously  it 
is  not  difficult  to  investigate  the  wave  resistance  corresponding  to  any 
curve  of  the  family  defined  by  Equation  [4e]  at  the  three  depths  of  immersion 
for  which  the  Integrals  have  been  tabulated. 

*There  it  •  slight  error  in  the  rteiattnco  curvet  R  of  Model  III  due  to  inaccuracy  in  computations, 
but  it  dost  not  invalidate  the  ooapariton. 
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Figure  20  -  Wave-Resistance  Coefficients  cw  =  ig  Referred  to  the  Wetted 

Area  S  for  the  Four  TMB  Models  I-IV,  f/L  =  0.125  (For  Comparison 
of  Range  Single  Curves  for  f/L  =  0.25  and  f/L  =0.5  are  Shown) 

To  check  the  order  of  magnitude  of  the  wave  resistance  and  to  en¬ 
able  a  comparison  with  experimental  data,  resistance  coefficients  c  of  the 

w 

four  TMB  models  I  to  IV  are  shown  in  Figures  20  to  22,  calculated  for  b/a  =  1  /~f 
and  C  =  1.07-  In  this  case  the  depth  of  immersion  ratios  f/L  correspond  to 
the  technically  more  familiar  f/D  ratios  as  follows: 


f/L 

0.125 

0.25 

0.5 

f/D 

0.875 

1-75 

3-5 

Assuming  a  rather  high  viscous-drag  coefficient  (cy  =  0.005),  the 
relative  importance  of  the  wave  resistance  at  various  depths  of  immersion 
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Figure  21  -  Wave-Resistance 
Coefficients  c„  =  — ,—IL. — 
w  p/2  VS 

Referred  to  the  Wetted  Area 
S  for  the  Four  TMB  Models 
I-IV,  f/L  =  0.25 


and  Froude  numbers  can  be  estimated  for  a  comparatively  wide  span  of  prismatic 
coefficient  0.71  ^  ^  0.59*  Attention  is  drawn  to  the  changes  in  the  mutual 

relations  between  the  curves  in  Figures  20  to, 22.  These  changes  are  dependent 
upon  f /L  and  upon  the  obvious  shift  of  the  peaks  towards  smaller  Froude  numbers 
as  compared  with  Figures  16  to  18,  because  of  the  factor  U2  in  the  denomina¬ 
tor  of  c„. 
w 

Considerations  of  wave  resistance  may  play  a  significant  role  when 
fixing  the  optimum  elongation  ratio  D/L  as  long  as  free-surface  conditions 
are  important.  Assuming  both  V  and  <t>  to  be  constant,  the  surface  S  and  there¬ 
fore  the  viscous  drag  vary  only  with  Vl7d  while  the  wave  resistance  varies 
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wit.u  (D/L)2  multiplied  by  a  complicated  function  rQ  of  P.  Restricting  P  to 
a  range  ~0.6  ^  P  ^  0.35.  rQ  is  monotonically  and,  on  the  average,  heavily 
decreasing  with  decreasing  P.  Thus  any  reduction  of  D/L  heavily  reduces  the 
wave  resistance. 

4.4.  LIFTING  DEPTH  OP  IMMERSION 

It  is  important  to  know  below  what  depths  of  immersion  f  the  wave 
resistance  can  be  neglected.  This  limit  can  be  established  from  such  cross 
curves  as  shown  on  Figure  11;  it  obviously  depends  upon: 

a.  The  Froude  number  P  or  yQ  =  1/2F2 

b.  The  L/D  ratio,  and 

c.  The  dimensionless  shape  of  the  body,  primarily  its  prismatic  coef¬ 
ficient  <(>  ,  especially  outside  of  the  large  hump. 

However,  some  additional  simple  reasoning  may  be  helpful  when  curves 
R  =  R(f/L)  are  not  available.  We  can  consider  the  wave  resistance  as  negli¬ 
gible  either  when 

a.  It  is  a  small  percentage  of  a  given  standard  resistance,  or 

b.  It  is  less  than  an  absolute  small  value  4R. 

Some  obvious  differences  in  results  due  to  the  different  approach 
have  sometimes  been  overlooked.  J 

a.  Assume  that  for  f  >  f  the  wave  resistance  becomes  less  than  a 

given  small  fraction  e  of  the  wave  resistance  RQ  at  a  standard  depth,  for 

instance  at  the  immersion  of  one  diameter;  f„  is  derived  from  a  ratio  of  the 

o 

resistances  in  question.  Comparing  bodies  of  equal  length,  fQ  depends  upon 
the  Froude  number  and  upon  the  dimensionless  shape  of  the  body,  but  only  very 
slightly  upon  the  elongation  ratio  D/L  ss  b/a,  since  the  latter  influences 
only  the  constant  ^jrC2pg  bVa,  which  drops  out  in  the  comparison. 

b.  Assume  that  the  limiting  depth  fQ  is  derived  from  the  condition 
that  the  wave  resistance  is  less  than  an  absolute  value  6R  independent  of 
the  standard  resistance  RQ.  Comparing  again  bodies  of  equal  length  fQ  now 
becomes  highly  sensitive  to  changes  in  D/L. 

A  rough  idea  of  the  necessary  limiting  depth  f  of  immersion  can 
be  obtained  from  the  decline  of  the  water  disturbance  with  increasing  depth 
in  a  plane  sinusoidal  wave;  this  estimate  normally  gives  exaggerated  values 

V 

Denoting  the  wave  amplitude  by  nm  and  the  amplitude  of  the  distur¬ 
bance  by  h  one  obtains 

-111 
k 


h  ■  \ e 
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putting  furtner 


x  =  Mf  =  2n?zL 


g 


h  =  h  e 
ci 


ik 


and  prescribing  h/hm>  for  Instance  assuming  h/hm  <  0.01  ,  one  obtains 


fo  >~0.75* 
or 

fQ/L  >  1.5?rF2 

This  estimate  Is  superficial  for  many  reasons: 

a.  The  resistance  depends  rather  on  the  square  of  the  generated  wave 
amplitudes , 

b.  The  actual  problem  is  three  dimensional,  and 

c.  The  body  shape  is  neglected. 

However,  it  shows  at  least  that  in  principle  the  limiting  depth  cannot  be 
expressed  as  a  fraction  of  the  dimensions  of  the  body  alone,  since  it  depends 
upon  the  length  of  the  free  wave  A  or  the  Proude  number  P. 

From  practical  considerations  matters  are  somewhat  different.  As 
mentioned  before,  at  very  high  Proude  numbers  the  ratio  of  wave  resistance 
to  frictional  drag  is  normally  very  small.  Thus  the  problem  of  finding  an 
accurate  value  of  the  limiting  depth  becomes  rather  unimportant  since  even 
grave  errors  in  computing  it  do  not  lead  to  appreciable  errors  in  the  total 
resistance. 


5.  BODIES  OF  REVOLUTION  OP  LEAST  WAVE  RESISTANCE 
5.1.  TWO-PARAMETER  FORMS 

In  an  earlier  paper5  endeavors  were  made  to  determine  distributions 
of  least  resistance  for  given  Proude  numbers.  The  results  varied  with  Froude 
numbers  and  depths  of  immersion,  which  is  quite  natural  in  the  light  of  such 
resistance  graphs  as  represented  by  Figures  6,  7  and  8. 

An  important  feature  is  the  peculiar  "swan  neck"  form  obtained  for 
higher  Froude  numbers -equal  to  and  above  P  =  0.35-  Because  of  the  limited 
accuracy  of  these  former  calculations  the  problem  has  been  reconsidered  here. 
The  present  investigation  supports  the  earlier  statements. 

The  formalism  needed  is  very  simple.  Some  controversy  arose  as  to 
how  far  the  application  of  exact  methods  of  the  calculus  of  variation  is 
consistent  when  dealing  with  surface  ships;5  the  results  obtained  did  not 
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lead  to  reasonable  ship's  forms.  However,  when  we  restrict  ourselves  to  fam¬ 
ilies  of  curves  expressed  by  polynomials  with  few  arbitrary  parameters,  we 
really  obtain  an  ordinary  minimum  problem  and  do  not  need  to  bother  about  the 
difficulties  connected  with  the  application  of  the  calculus  of  variations. 
Take  for  instance  the  family  (basic  form) 

u  -  l  -  -  a -  <‘)  -  a -  (•)  (34) 

with  two  arbitrary  parameters.  The  wave  resistance  R  is  given  as  a  second 
degree  function  in  a  and  a  . 

2  4 

R  =  4B  a2  +  4B  a2  +  8b  a  a  +  24B  a  +  24B  a  +  Bn  [35] 

22  2  44  4  24  24  22  44  0 


where 


B  =  ft  -  6ft  +  9ft. 

22  11  15  55 

B  =  4731  -  1  277?  +  9ft 

44  S3  35  55 

B  =  2ft  -  3ft  +  9ft  -  6ft 

24  13  Is  55  3! 


B  =  ft  -  3ft 

2  15  55 

B  =2 ft  -  3ft 

4  35  5! 


differentiating  R  partially  with  respect  to  a^  and  a^,  one  obtains  the  min¬ 
imum  conditions 

P-  =  B  a  +  B  a  +3B  =  0 

oa  22  2  24  4  2 

2  [36] 

P-  =  B  a  +  B  a  +3B  =  0 

o a  24  2  44  4  4 


whence 


3[B  B  -  B  B  ] 

a  =  2  44  4  24 

2  B  B  -  B2 

22  44  24 

3  [B  B  -  B  B  ] 
a  =  4  22  2  24 

4  B  B  -  B2 

22  44  24 


These  equations  lead  to  results  which  are  not  applicable  to  practice 
when  yQ  =  1  and  of  restricted  interest  when  yQ  ■  2  (f/L  is  assumed  equal  to 
0.125). 
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Figure  23  -  Doublet  Distribution  for  Least  Resistance, 

Two-Parameter  Forms,  F  =  0.408,  Yq=  3 

When  y Q  -  3  the  distributions  shown  on  Figure  23  are  obtained.  We 
note  again  the  difference  in  the  shapes  when  f/L  =  0.125  and  f/L  =  0.25.  Ex- 
tending  the  calculations  to  =  4  and  y0  =  5>  curves  of  more  and  more  "rea¬ 
sonable"  character  are  obtained  as  shown  in  Figures  24  and  25. 

The  apparent  failure  of  the  theory  to  yield  useful  results  in  some 
cases,  is  often  due  to  lack  of  suitable  conditions  imposed.  There  is  no 
reason,  for  instance,  to  expect  a  solution  which  leads  to  a  "normal"  prismatic 
coefficient  if  no  restrictions  as  to  this  coefficient  are  made.  On  the  con¬ 
trary,  it  is  rather  fortunate  that  one  obtains  results  which  meet  other  re¬ 
quirements  of  practice  (i.e.,  are  "reasonable"),  without  this  restriction  in 
certain  ranges  of  Froude  numbers. 

5.2.  ISOPERIMETRIC  PROBLEMS,  ONE-PARAMETER  FORMS 

Introducing  a  condition  0  =  const  we  obtain  an  isoperimetric  pro¬ 
blem.  Then  Equation  [34]  retains  only  one  arbitrary  parameter.  This  can  be 


I  1  I  I  I  I  i 

f/L  «Ol25t  i»  «  <-  3.I97|*  ♦  6.667 14  -  4.470#*,  4  *  0.628,  t  »  6.55 
vJ/L »  0.25  ,  9  *  I  “4 .617 1*  ♦  *0.915  #4  -  7.298  #*,  4  »  0.601,  t  *  8.36" 

.  X  I  I  I  I  I  I  I 
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interpreted,  for  example,  as  Taylor's  tangent  value  t.  The  resulting  equation 
is  of  the  type 

*(0  =  >70U)  -|t'  [*2  -  ^  f4  +  ^#]  [38] 

Here  »70(()  is  a  given  polynomial  complying  with  the  conditidn  <j>  -  const;  its 
tangent  value  tQ  may  be  chosen  in  such  a  way  that  the  equation  rjQ  is  as  simple 
as  possible.  The  function 

-|[f2  <4  +  j  <‘]  =  [38a] 

has  the  properties: 

34  nO) 

1  •  “a?  =  "1 

2.  4  rj  (o)  =  4  rj{  1  )  =  0 

2  2 

3.  f4  nUJdf  -  0 
Jo  2 

t'  is  the  variable  tangent  parameter,  the  resulting  t  of  the  Equation  [38] 
being  obviously  t  =  t'  +  t  . 

Assuming  0  =  2/3 ,  rjQ  «  1  -  t 2 


=  -2f  -  |t'  [(  -  +  7(=] 

one  obtains 

R  -  t,2Aa  +  3t'Ai  +  4  7^  [39] 


aR  _  9 
at>  "  j 


t'A 

2 


+  3A  ■  0 


[39a] 


with 


15 


280 
—  Kt 


3b 


rr\ 


11 


TTO 


+  77H 

13  15 


hence 


[39b] 


Another  isoperiraetric  problem  is  given  by  t  «  const  and  variable. 
Although  this  problem  looks  somewhat  artificial  since  there  are  no  technical 
reasons  to  keep  the  tangent  of  the  sectional -area  curve  rigidly  fixed  the 
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Figure  26  -  Doublet  Distribution  for  Least  Resistance, 
One -Parameter  Forms,  F  «=  0.408,  7Q  «  3 


Figure  27  -  Doublet  Distribution  for  Least  Resistance, 
One-Parameter  Forms,  F  =  0,408,  YQ  =  3 
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results  are  Interesting.  Figures  26  and  27  show  that  assuming  rather  dlffer- 
ent  t  values  optimum  ship  lines  with  a  similar  trend  may  be  obtained. 

We  notice  that  the  optimum  area  coefficient  *  for  a  medium  depth  of 
immersion  f/L  *  0.25  is  much  higher  than  for  a  slight  immersion  f/L  *  0.125. 
This*  might  have  been  inferred  from  the  shift  of  the  resistance  curves  follow¬ 
ing  Figures  6,  7  and  8. 

The  necessary  formalism  is  again  very  simple:  assuming  as  before 
a  curve  tjq  with  the  fixed  t  =  tQ  value,  for  example  as  before,  17 Q  =  1  -  f2, 
tQ  *  2,  and  denoting  *  =  *0  +  *' » 


n  •  1  -  *2  +  13.1 25<*' U2  -  2£4  +  ?) 
&xnU)  =  13.125U2  -  2*4  +  {•) 

complies  with 

1.  A  ij(0)  =  4^(1 )  =  0 

2.  \x  A  ijU)d*  =  1 
dA  if(l  ) 

3.  -^—-o 

from 

-  -2*  +  26.25*'  l*  -  4*3  +  3*5] 

we  obtain 

2  A' 

*'  "  205  AT 
2 


with 

A'  =  hi  +  l6to  +  9771  -87M  +  67T1  -  2474 

2  11  33  55  13  i 5  35 


[40] 

[40a] 


[41] 


A' 

x 
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6.  RESISTANCE  CURVES  OF  THE  FAMILY  (2,  4,  6;  * ;  t) 

A  systematic  survey  of  resistance  properties  of  ship  forms  can  be 
obtained  by  a  different  approach,  i.e. ,  varying  the  parameters  of  a  given 
family  of  ship  lines  and  plotting  the  corresponding  resistance  curves.  Re¬ 
stricting  ourselves  to  a  two-parameter  equation 

(2,  4,  6;*;  t)  -  1  -  £  a^* 

2,4,6 


[34] 
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Equation  [35]  can  be  used  for  calculating  the  resistance,  or  still  simpler, 

R  *  4[a2  it?  +  4a  *Tt\  +  9a2W|  +  4aa  «I  +  6a  a  +  12a  a  7T)  ] 

2  11  4  33  C  35  2  4  13  2  1  13  4  <  33 

[42] 

The  parameters  a  ,  a^  and  a#  are  connected  with  the  basic  form  co¬ 
efficients^  and  t  by  the  equations 


a  *  1  -  a  -  a 

*  2  4 

Table  3  contains  wave -resistance  coefficients  rQ  for  t  »  0,  1,  2,  3  and 
0.68  >  ^  >  O.56  with  an  interval  of  A+  ■  0.02  within  a  range  of  Froude  num¬ 
bers  1  >  F  >0.25  (for  t  *  0  additionally  *  0.50,  O.52,  0.54)  at  a  depth  of 
immersion  ratio  f/L  ■  0.125. 

The  corresponding  curves  spaced  A+  »  0.04  are  shown  on  Figures  28 
to  35  grouped  following  t  and  4.  The  main  purpose  of  these  plots  is  to  dem¬ 
onstrate  -the  dependence  of  the  wave  resistance  upon  t  for  0  =  const;  it  is 
interesting  to  note  that  the  peak  values  (cf  page  24)  differ  as  much  as  -15 
percent  for  t  «  0  and  t  «  3,  in  close  agreement  with  results  known  from  stu¬ 
dies  of  surface  ships  and  the  tendency  exposed  by  the  minimum  calculations. 

One  should,  however,  remember  that  theory  tends  to  overestimate  the  favorable 
interference  effects  and  that  viscosity  precludes  the  realization  of  excessive 
angles  of  run.  On  the  other  hand,  for  very  high  Froude  numbers  the  relative 
importance  of  asymmetry  decreases,  so  that  forms  with  steep  slopes  at  the  bow 
and  moderate  slopes  at  the  stem  may  be  advantageous. 

SUMMARY 

Using  Havelock's  basic  work  and  some  former  investigations  by  the 
present  author,  a  systematic  synopsis  is  made  on  the  wave  resistance  of  bodies 
of  revolution.  Tables  evaluated  by  the  Bureau  of  Standards  and  graphs  are 
given  which  allow  the  investigator  to  estimate  immediately  the  wave  resistance 
of  a  wide  class  of  bodies  of  revolution  defined  by  doublet  distributions  along 
the  axis  expressed  by  polynomials. 

Some  discussions  refer  to  the  relations  between  this  distribution 
fx*(0  and  the  sectional  area  of  the  body  A#U)*  For  "normal"  shapes  of  dis¬ 
tribution  the  usual  assumption  is  made  that  there  is  affinity  between  /**({) 
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and  A#({).  In  extreme  cases  the  shape  of  the  body  can  be  calculated  by 
methods  due  to  Landweber  and  Amtsberg;  no  corrections,  however,  are  given  for 
the  Influence  of  the  free  surface  on  the  shape. 

Within  the  first-order  theory  the  resistance  can  be  split  up  into 
a  main  part  due  to  a  symmetric  distribution  with  respect  to  the  midship  sec¬ 
tion  and  a  part  due  to  asymmetry,  which  can  be  investigated  Independently. 
Large  amounts  of  asymmetry  can  influence  the  resistance  detrimentally  in 
some  ranges  of  the  Proude  number. 

The  investigation  of  the  resistance  as  a  function  of  the  body  form 
leads  to  conclusions  which  sometimes  are  contrary  to  those  derived  for  sur¬ 
face  ships.  The  choice  of  appropriate  prismatic  coefficients  varies  deci¬ 
sively  with  the  range  of  the  Prouue  number,  as  is  clearly  illustrated  by 
the  numerous  graphs.  The  same  applies  to  the  influence  of  the  tangent  value 
t.  Ceteris  paribus  the  resistance  is  approximately  proportional  to  the 
square  of  the  midship  section. 

The  dependence  of  the  resistance  upon  the  depth  of  immersion  is  ln-^ 
vestlgated;  this  dependence  is  best  explained  by  the  ratio  f A,  where  A  = 
is  the  length  of  the  free  wave.  Thus  for  conmon  prismatic  coefficients  the 
wave  resistance  decreases  rapidly  with  increasing  f  except  in  the  range  of 
high  Proude  numbers  (large  A  values).  In  the  range  of  high  P  the  calculation 
of  forms  (distributions)  of  least  resistance  leads  sometimes  to  results  bare 
of  practical  applicability;  by  Introducing  suitable  restrictions  such  diffi¬ 
culties  are  avoided.  Wiese  investigations  show  important  peculiarities  of 
the  distributions. 

A  set  of  resistance  diagrams  calculated  for  the  family  (2,4, 6;^, *t) 
gives  a  survey  of  the  resistance  properties  of  a  class  of  normal  bodies. 
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Figure  28  -  Wave-Resistance  Coefficients  rQ  =  ij<rcrpg''b'r7a» 
of  Synnnetrical  Bodies  Belonging  to  the  Class 
(2 ,4, 6;0; t )  «  1 -a2$2-a4$4-a#$4  Over  YQ  -  1/2P2  and  F  »  V/lfiL 


Figure  29  -  Wave -Resistance  Coefficients  rQ  ■  4^c'*pg  ^  /a“» 
of  Symmetrical  Bodies  Belonging  to  the  Class 
(2,4,6;$;t)  *  1-a  *2-a  f*-a  l*  Over  Y  «  1/2F2  and  F  ■  U/fSIT 

2  4“  3 

t  ■»  1 


Figure  32  -  Wave -Resistance  Coefficients  r  =  ^  R  t, 

of  Symmetrical  Bodies  Belonging  to  the  Cl£ss8  b 
(2,4,6;*;t)  =  l-a^2-a^4-a^6  Over  yQ  =  1/2P2  and  P  «  u/j^L 

0  *  O.56 


Figure  33  -  Wave -Resistance  Coefficients  rQ  ■  jr — y  R — T7“* 
of  Symmetrical  Bodies  Belonging  to  the  Class8  *  * 
(2,4,6;^;t)  =  1  "a2*2"a4*4-a6*#  Over  yQ  =  1/2F2  and  F-«  U /VgL 

<t>  =  0.60 
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APPENDIX  I 

APPROXIMATE  CALCULATION  OP  THE  SURFACE  S  OF 
A  CLASS  OF  ELONGATED  BODIES  OF  REVOLUTION 

From  Guidin' s  rule 

S  =  lit  |  yds  *  2ir  £“ y  /l  +  (y1)2  dT 
S  s  2n  J+*ydx  +  it  J**y  y'2  dx  ■  wAn  +  n J  *  y  y,a  dx  [29] 


the  main  part  of  the  surface  is  given  by  n  times  area  of  the  meridian  section 
Am  plus  a  correction  term  neglecting  higher  order  terms. 


With  y  «  bH 


y 


_b  dH 
a  Ff 


the  correction  term  becomes 


IT 


fyy 


|2 


dx  «=  n-ab 


d*  -  »rab  ^1  .  [30] 
a* 


i.e.,  the  correction  term  is  equal  to  the  area  of  an  ellipse  with  the  axes 
a,  b  multiplied  by  the  square  of  the  elongation  ratio  and  a  numerical  value 
I  dependent  upon  the  equation  of  the  curve.  To  get  an  idea,  with  obvious 
denotations, 


I2  -  16/15 

H  -  1  - 

I4  -  128/77 

Mr 

1 

n 

X 

T  2n® 

H  8  1  * 

Tn  "  (2n-1  T  (3n-fJ 

The  next  term  in  the  expansion  of  S 

'  f  C  y  y'4  dx  -  ^  J*1  HH'4d*  [31  ] 

fit 

with  H'  -  dH/d*  is  obviously  of  the  order  b4/a4.  However,  taking  H  »  1  -  $* 
the  factor 

J>‘  «  -  K  ■ 

grows  with  n3  when  n  is  large.  Provided  b/a  is  not  too  small,  say  ~1/7,  the 
error  in  neglecting  all  terms  except  the  fir.st  (Equation  [29])  is  only  per¬ 
missible  as  long  as  n  <  -5. 
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Using  Equation  [30],  various  changes  In  S  can  be  easily  estimated 
within  the  range  of  validity  of  the  formula.  Por  Instance,  the  Influence  of 
an  asymmetric  term  can  be  discussed  as  follows: 

Ht  *  Hs  +  "a 

h  '  £‘  <»,  +  H a)(Hi  +  '  I,  +  £'  HsH^d<  +  £'  d{ 

[32] 

'Where  I„  refers  to  the  even  part  following  [30].  When  the  meridian 
curve  has  vertical  tangents  at  the  bow  and  stern  (or  bow  or  stern)  the  pre¬ 
ceding  reasoning  can  be  applied  In  principle  to  a  range  1  -  >{>-(1  -  «A.), 

and  the  remainder  Is  calculated  as  the  surface  area  of  a  segment  of  the 
sphere  generated  by  the  radius  of  curvature  at  the  nose  or  stern.  Such  an 
approach  is,  however,  only  useful  when  the  Integrals  Involved  are  of  simple 
type. 


m 
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APPENDIX  II 

EVALUATION  OP  THE  AUXILIARIES  INTEGRALS* 
The  Integral  to  be  computed  is  given  by: 

r-  .*£?  (y/yJ2 

W  _  t  _  r*  I  .  O  _  u  /  ..\u  /..I 


I  =  c 


. M.fyllMyJdr 

(y/y0)2-i  1  J 


The  functions  e  and  M^yjM.jfy)  are  well-behaved  in  the  entire  interval  of 

j x.  i_  i a  « xi  _  1  _  ft a  _ i  » a  /  ■  /O ..  <  _  _  _  .  


integration.  However,  the  algebraic  function  (y/y#)/  V(Y/yt)* -  1  causes  some 
difficulty. at  the  lower  limit  of  integration, ' i.e. ,  at  y  *  yQ=  In  the  neigh¬ 
borhood  of  y  *  yQ,  the  contribution  of  I  is  far  from  negligible  and  therefore 
an  investigation  was  carried  out  to  determine  the  asymptotic  behavior  of  the 
integral  as  a  function  of  the  upper  limit.  Specifically,  the  following  func¬ 
tion  was  examined: 

+ 


>  -  r 

Jr . 


M(y)f(y)dy  «  >  0 


where 


M(y)  -  e  V MjtrIMjlr) 

Mr.)2 

w_.\  -  _  o 


It  was  found  that: 

I(«)  -  e~*r°  M^(y )Mj ( y)  yQ  Yu  |l  +  +  0(«2)| 

This  asynqjtotlc  expression  was  used  to  determine  the  interval  of 
integration,  Ay  far  a  numerical  integration.  This  Interval  was  too  small  to 
be  practicable,  even  allowing  for  subsequent  changes  in  Ay. 

A  new  approach  to  the  problem  was  sought  in  u  suitable  transforma¬ 
tion.  The  following  trs ns format ion  very  quickly  presented  itself: 


y  ■  Z2  +  yQ 
dy  =  2Z  dZ 

The  original  Integral  was  transformed  as  given  by: 


*0?  J.  BIub,  fatiooal  Buratu  of  Standards 


ft- 
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£ 


1  -  ~  f  e~  ^  M1  (y)Mj(  y)  (Z2  +  y )dZ 

ro  *  Vz *+2*.  1  J  0 


In  this  form,  the  integrand  behaves  properly,  (there  is  no  longer 
a  singularity  at  y  *  yQ)  and  the  integral  converges  rapidly. 

The  Integral  was  actually  confuted  by  using  the  fora  in  [5].  The 
numerical  integration  was  performed  once  using  Simpson's  rule  and  a  second 
time  using  the  trapezoidal  rule— for  checking  purposes.  The  Interval  dZ  was 
taken  as  0.1  and  the  range  extended  from  0.0  to  approximately  3.5*  The  M 
functions  vere  computed  from  previous  tables  by  using  4-point  Lagranglan  In¬ 
terpolation.  The  exponential  function  was  computed  from  tables  and  the  use 
of  the  approximation  e”x  «  1  -  x  +  xa/2,  x  <0.01.  The  algebraic  function 
in  the  Integrand  was  computed  in  straight  form  and  fashion.  All  of  this  work 
was  done  on  the  IBM  electronic  calculator  (type  604)  and  the  auxiliary  IBM 
punch  card  equipment.  All  the  IBM  operations  were  checked—  Independently 
wherever  possible. 


Auxiliary  Integrals  for  Various  Proude  Numbers 
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.000485850 

0.000601895 

0.010330? 

.00611066 

.00376711 

.00329053 


.00470793 

.00799232 

0.0131527 


0.00150426 

.000527452 

.OO13C715 

.00384335 


.0081  3606 
.0141853 

0.0219910 


0.000423770 
. 000782298 
.00255238 
.00573401 


.0103272 

.0163319 

0.023748? 


0.000271083 
. 00081 6457 
.00228089 
.00466439 


.00796694 

.01  21886 
0.017329? 


0.000301165 
.  00047881 7 
.00114028 
I  .00228557 


.00391467 

.00602758 

0.00862431 


0.000553898 

000428905 

.000500888 

.000769848 


.00123578 
.  001 89870 
0.00275892 


0.000882234 

.000746938 

.000667035 

.000642525 


. 000673406 
.  000759681 
0.000901347 


0.00102788 
.00104870 
.001 08303 
.00113087 


.00119221 
.00126705 
0.001  35540 


0.00884174 

.00468247 

.00239918 

.00199186 

0.00185989 
. 000284792 
. 000466200 
.00240411 

.00346053 
.00680518 
o.oi  20258 

.00609853 

.0115494 

0.0167569 

o.ooi  78332 

.000947282 

.00152280 

.00350987 


.00690849 

-0117187 

0.0179404 


O.OO181881 

.001051  65 
.001 20355 

.00227450 

0.00246152 

.00156282 

.00114793 

.00121685 

0.00344350 
.  00264391 
.  002C41 29 
.ooi  ^3565 

0.003843?? 

.0034?7?8 

.00306673 

.00276158 

0.00331561 

.00330874 

.00331538 

.0033355? 

.00426453 
.0071 7361 
0.0110018 

.00176959 

.00280615 

0.0043265? 

.00142698 
.  001 41 529 

0.001 6005P 

.00251182 

.00231745 

0  00217848 

.00336916 

.00341631 

0.00347697 
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